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Solutions of Laplace’s equation with simple boundary conditions, and their 
applications for capacitors with multiple symmetries 

Mayckol Morales* Rodolfo A. Diazj William J. Herrera! 

Departamento de FIsica. Universidad Nacional de Colombia. Bogota, Colombia. 


Abstract 

We find solutions of Laplace’s equation with specific boundary conditions (in which such solutions take either the value 
zero or unity in each surface) using a generic curvilinear system of coordinates. Such purely geometrical solutions (that we 
shall call Basic Harmonic Functions BHF’s) are utilized to obtain a more general class of solutions for Laplace’s equation, 
in which the functions take arbitrary constant values on the boundaries. On the other hand, the BHF’s are also used to 
obtain the capacitance of many electrostatic configurations of conductors. This method of finding solutions of Laplace’s 
equation and capacitances with multiple symmetries is particularly simple, owing to the fact that the method of separation 
of variables becomes much simpler under the boundary conditions that lead to the BHF’s. Examples of application in 
complex symmetries are given. Then, configurations of succesive embedding of conductors are also examined. In addition, 
expressions for electric fields between two conductors and charge densities on their surfaces are obtained in terms of gen¬ 
eralized curvilinear coordinates. It worths remarking that it is plausible to extrapolate the present method to other linear 
homogeneous differential equations. 

Keywords: Laplace’s equation, curvilinear coordinates, capacitance, harmonic functions, separation of variables. 
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c/> Introduction 

Solutions of Laplace’s equation, usually called Harmonic Functions (HF’s) are very important in many branches of Physics and 
: '"’T ingi peering such as electrostatics, gravitation, hydrodynamics and thermodynamics[I] [2]. Therefore, considerable effort has 
been done in solving Laplace’s equation pF|-p?5] in a variety of geometries and boundary conditions. On the other hand, the 
capacitance of electrostatic conductors is a very important quantity since capacitors are present in many electric and electronic 
devices m-m- Thus, many studies of its general properties m-m, and calculations of coefficients of capacitance for many 
QQ geometric configurations [T5]-p?8] have been carried out. 

t—H The first goal in this paper is to show a simple method to solve Laplace’s equation for configurations of volumes with certain 

0“N symmetries, when the HF’s acquire values of either zero or unity in each surface that provides the boundary conditions. The 
CD solutions of Laplace’s equation under these specific boundary conditions will be called Basic Harmonic Functions (BHF’s). 
Cf) We shall see that the form of these boundary conditions leads to very easy solutions via separation of variables. Besides, the 
l w ) solutions will be written in terms of a generic system of curvilinear coordinates that can be adjusted for many symmetries. It 
is remarkable that the BHF’s depends exclusively on the geometry, and not on the physical problem involved. 

. . The determination of the BHF’s under a given geometry will be applied in two scenarios: (a) The generation of a more 
general class of solutions in which the HF’s take arbitrary constant values on the surfaces that determine the volume, (b) The 
calculation of the coefficients of capacitance when the surfaces that form the volume are covered with electrostatic conductors. 

< i Once again, both types of solutions are given in a generic form in terms of appropriate generalized systems of coordinates. 

03 The paper is distributed as follows: In section [1] we solve Laplace’s equation in generalized geometrical configurations to 
obtain the BHF’s, by using a generic system of curvilinear coordinates. As a matter of illustration of the method, we obtain 
the BHF’s for the cases of two concentric spheres and two concentric cylinders. Then, from the BHF’s obtained in Sec. |T]and 
appealing to the linearity of Laplace’s equation, we construct in Sec. [2] a more general class of solutions in which HF’s acquire 
arbitrary constant values on the surfaces. On the other hand, in section [3] we use the general expressions for the BHF’s in 
generalized curvilinear coordinates, in order to obtain the corresponding generalized formulas for the coefficients of capacitance. 
The capacitance coefficients for the cases of two concentric spheres and two concentric cylinders are obtained for illustration. 
Then, section [4] shows calculations of BHF’s and capacitances for more complex geometries in which our formulation acquires 
all its power. Moreover, in appendix [A] we obtain the electric field between two conductors as well as charge densities on their 
surfaces in terms of curvilinear coordinates, while in appendix [Bj formulas for the coefficients of capacitance are extended to 
the case in which we have a configuration of succesively embedded conductors. Appendix [C] provides some relations between 
curvilinear coordinates and geometrical factors. Finally, section [5] contains our conclusions. 
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1 Laplace’s equation and systems of orthogonal curvilinear coordinates 


Laplace’s equation given by 


V 2 / = 0 


(1) 


has solutions that depend on the boundary conditions. It is desirable to use a system of coordinates adjusted to the symmetry of 
the problem. Let us consider a generic orthogonal coordinate system (it, v, uj), with its corresponding scale factors ( h u , h v , h w ). 
Furthermore, let us assume that the boundary consists of two closed surfaces S\ and S 2 in which the coordinate u takes 
constant values u\ and 112 on S\ and S 2 respectively. We shall suppose that one of the surfaces contains the other, and we 
define as V the volume delimited between both closed surfaces. We intend to find a solution fj(u,v,w ) of Laplace’s equation 
in V, with the following boundary conditions 


fj (uj,v,w) = 1, fj(ui,v,w) = 0 ; with i,j = 1,2 and i j 


( 2 ) 


where i, j = 1,2 are indices that label the surfaces that determine the boundary. Further, Ui refers to the constant value of the 
coordinate u on the surface Si. As we already mentioned, we shall call Basic Harmonic Functions (BHF’s) to the solutions of 
Laplace’s equation with the boundary conditions defined in Eq. ©• In terms of the ( u , v, w) coordinates, Laplace’s equation 
for fj becomes [T] 


d_ f h v h w dfj \ 9 f hwhu dfj 
du \ h u du) dv \ h v dv 


d f h u h v dfj \ 
dw \ h w dw J 


= 0 


( 3 ) 


Assuming the following ansatz of separation 


fj (u, v, w) = Uj(u)Qj(v, w) ; j = 1,2 (4) 

The boundary conditions Eq. m yields 

fj ( Uj,v,w ) = Uj ( Uj) &j(v,w) = 1 ; j = 1,2 (5) 

fj(ui,v,w) = Uj(ui)Qj(v,w) = 0 ; with i,j = 1,2 and i j (6) 


From now on, we shall assume that i,j = 1,2 and i ^ j unless otherwise stated. Since the coordinates v and w take any 
value on the surfaces, it is necessary that Uj ( m ) = 0, in order to satisfy the boundary condition ([6j and to obtain a non-trivial 
solution. On the other hand, the boundary condition (0 says that 


0j (v,w) 


1 

Uj (uj) 


= Ai 


( 7 ) 


where Aj are two constants that only depend on j. Thus, the solution fj(u,v,w) in Eq. 0 can then be rewritten as 


fj(u, v, w ) = AjUj (u) = fj (u ) 


( 8 ) 


therefore the solution only depends on the coordinate u. The boundary conditions 0 become 


from Eq. @ we get d v fj = d w fj 


fj ( m) = 0 ; fj (■ Uj ) = 1 
0 from which Laplace’s equation 0 reads 


d_ / h v h w dfj (u) \ 
du \ h u du J 
h v h w dfj (u) 
h u du 


0 => 
kj (v, w) 


( 9 ) 


( 10 ) 

( 11 ) 


note that there is not dependence on the coordinate u on the right hand side of Eq. HD- Consequently, the u— dependence 
given by d u fj ( u ) on the left hand side of Eq. CD must be cancelled by the remaining factor h v h w /h u . It suggests that the 
dependence with u on the term h v h w /h u should be factorized. Therefore, it is reasonable to assume that 


h V hy 


= G(u)H (v,w) 


( 12 ) 


Note that the validity of this ansatz depends only on the coordinate systeix0. Thus, neither G ( u ) nor H (v, w ) could depend 
on j . Applying such an ansatz on Eq. (1101) we obtain 


In practice, the ansatz of separation 02} holds in many systems of curvilinear coordinates. 
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H (v, w) 


du 


G(u) 


dfj (u) 

du 


= 0 => 


d 

du 


G{u) 


dfj {u) 


du 


= 0 




where Bj is another couple of constants for each Sj. We can solve for fj (u) as follows 


dfj (u) 
fj (u) 


Bj 

G(u) 


du => 


B j Z{u) + C j 


fj ( u ) ~ Bj f G{yL) + G i 

dZ (u) _ 1 

’ du G (u) 


(13) 

(14) 


(15) 


once again Cj is another couple of constants associated with each surface. Combining (fl5l) with the conditions (0 we have 



fj ( Ui ) = 0 = BjZ ( Ui ) + Cj 

(16) 


fj ( u j ) = 1 = BjZ ( Uj ) + Cj 

(17) 

From (flDl) it is obtained that 

Cj = —BjZ (uj) 

(18) 

Substituting (fl8l) in (flTl) we obtain 

1 = BjZ ( uj) - BjZ {m) 

(19) 


0 1 

J [Z ( Uj ) - Z ( Ui )\ 

(20) 


and substituting (fTDl) and EDI) in (fl5l) yields 


f ( v = Z(u) _ Z(u) _ Z {uj) 

Jl[ J [2W^W] [Z{u 5 )-Z{ui)] [Z( Uj )-Z( Ui )\ 

Summarizing, let us assume a geometric configuration with two closed surfaces Si and S 2 in which one of the surfaces 
contains the other. Both surfaces determine a volume V. Let us also assume that there is a system of coordinates (u,v,w), 
with scale factors (h u , h v , h w ) such that the coordinate u is constant in each surface S 1 and S 2 (both constants ui and should 
be different if we want a non-trivial solution). Under such hypotheses, the solutions of Laplace’s equation within the volume 
V and that satisfies the boundary conditions 


fj (uj,v,w) = 1 , fj(ui,v,w) = 0 ; with i,j = 1,2 and i yt j 


( 21 ) 


possess the following features: These solutions f j (u) only depends on the variable u and can be constructed with the following 
set of expressions 


fj (u) 


h v h 


W 


h u 


Z (u) - Z (■ m ) 
i z ( Uj ) - Z(ui)\ ’ 

G (u) H (v, w) ; 


with i,j = 1,2 and i yf j 

dZ (it) _ 1 

du G (u) 


( 22 ) 

(23) 


Note that 

fi (u) + /2 (u) = 1 (24) 

in the region V in which this Laplace’s equation is satisfied. Such a property is a particular case of a more general context [3j. 
From the previous developments, it can be seen that the functions G ( u ), H (v, w ) and Z (u) are not unique. We can redefine 
such functions as follows 

G (u) = kG (u) , H (v,w) =—H (v,w) , Z (u) = — Z (u) + fi (25) 

K K 

where k and /3 are arbitrary constants. Nevertheless, it is not a problem since these functions are not physical observables. 
Moreover, it is easy to check that the BHF’s given by Eq. f)22D are invariant under the “gauge transformations” defined by 
Eqs. E5D . Finally, the boundary conditions along with the uniqueness theorems, guarantee that each solution fj (u) is unique 
within the volume V. 

It worths remarking that the formulas obtained in this section are valid regardless which surface contains the other. Further, 
perhaps the most outstanding property of the BHF’s obtained here, is that they depend only on the geometry. The functions 
fj ( u ) are dimensionless and given a geometry, they are unique for each j. 
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1.1 BHF’s for a particular case for the scale factors 

Let us examine the particular case in which 

h v h w ^ 
hu 

in that case from Eq. (1DP1) the simplest solution for the functions G ( u), H (v , w) and Z (u ) reads 


(26) 


G (u) = H (v, w) = 1 ; Z (u) = u (27) 

Substituting (1PT1) in ED) the solution for fj ( u ) simplifies considerably 

U — 'LL ' 

fj (it) =-— ; with i,j = 1,2 and i ± j (28) 

'LL j 'LL'i 

The assumption (1261) is fulfilled when we have a system of coordinates with cylindrical symmetry that is obtained with a 
conformal transformation from the cartesian system of coordinates projected onto the XT—plane. In that case we usually 
obtain h u = h v and h w = 1 (of course the roles of v and w could be interchanged). Note however that a cylindrical symmetry 
coming from a conformal transformation is only a sufficient condition but not necessary. Equation (1261) is the only condition 
we require for the validity of Eq. ( 1251 ) . 


1.2 BHF’s for the case of two concentric spheres 

In order to illustrate the method, let us take a geometric configuration consisting of two concentric spheres of radius a and 
b with a > b. We recall that by definition u is the coordinate such that zz, keeps constant for all points of the surface Si for 
i = 1,2. Thus, the spherical coordinates are the appropriate ones, and u = r is the coordinate that is kept constant for each 
of the spherical surfaces. Hence, for this configuration we have 


from the first of Eqs. 


so a possible solution itQ 
from the second of Eqs. 


and substituting in Eq. E?j) we obtain 


u = r, v = 6, w = p and h r = 1, hg = r, h v = rsinf? 
we can find a solution for G ( u ) = G (r) and for H (v, w) = H (9, ip) 

= G (r) H (9,ip) => r 2 sin 9 = G (r) H (6, <p) 
h r 

G (r) = —r 2 ; H (9, ip) = H ( 9 ) = — sind 
we can find a solution for Z (u) = Z {r ) 


zw = -/ 


dr 


(29) 

(30) 

(31) 

(32) 


fj 0) = 


for rj = r 2 = a and = r\ — b with a > b Eq. 


i,j = 1,2 and i ^ j 


becomes 


f 2 M = 


h (r) = 


r i 


1 _ 1 
r b 


(i-i) (*-*) 

\^2 ri j 


ab 

b — a 


1 1 
r b 


(33) 


(34) 


since V 2 (1/r) = 0 for r ^ 0, it is easy to check that fa ( r ) is a solution for Laplace’s equation in the region b < r < a that 
satisfies the boundary conditions ([2]) that in our case read 

fi (r 2 ) = fi («) = 1 and f 2 (rq) = / 2 ( b) = 0 


^Note that G (r) = r 2 and H ( 6 ) = sin 0 is also a solution. Both solutions are connected by setting « = — 1 and j3 = 0 in Eqs. 1251 . We prefer 
solution ED to obtain a positive solution for Z ( u ). 
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from Eq. (1331) we can also obtain /i (r) 


h (7) 


h (r) 


1 1 


r 



ri 



(i 


1 

a 




(35) 


which satifies Laplace’s equation in the region b < r < a with boundary conditions / 1 (a) = 0 and fi (b) = 1. Observe that 
fi (r) + / 2 (r) = 1, as it must be. 


1.3 BHF’s for two concentric cylinders 

Let us assume two very long concentric cylinders of radii a and b with a > b. This geometry can be adapted to the polar 
cylindrical coordinate system (r,ip,z). However, we shall use the (equivalent) conformal polar cylindrical coordinate system 
(p, p, z ). Such a system is defined as 


x = e p cos ip ; y = e p sin tp ; z = z (36) 

h p = r ; hip = r ; h z = 1 (37) 


the advantage of using the coordinates (p, ip, z) [instead of the usual cylindrical polar coordinates] is that the scale factors (1371) 
satisfy the condition (1261) from which it is immediate the form of the BHF’s from Eq. (1751) 


f.( ) = p ~ Pi = ln 
3 P Pj - Pi In (rj/n) 


r\ = b ; r 2 = a with i,j = 1,2 and i/j 


(38) 


2 The case in which the surfaces are “equipotential” 

Returning to the general case, it could be realized that owing to the particular form of the boundary conditions (]7|). the 
solutions previously generated are remarkably simple and also purely geometrical as discussed above. We can use the BHF’s 
fj ( u) generated in this way, in order to obtain a more general class of solutions 0 ( u ) of Laplace’s equation in the same volume 
V, in which the HF 0 (u) acquires constant but arbitrary values 0 1 and 0 2 on the surfaces Si and S 2 respectively. For such a 
scenario, the solution 0 ( u ) is given by 

0(«) = (j)ifi (u) + <£ 2/2 ( u ) = 02 + fl {u) [01-02] (39) 

where we have used Eq. (1751) and the functions fj (u) are given by Eqs. (1771 1751) in generalized orthogonal curvilinear 
coordinates. We can see that Eq. (1571) is the solution to this “equipotential problem” by observing that if we apply V 2 on both 
sides of such an equation, we find that 0 (it) obeys Laplace’s equation. Moreover, by taking a point on the surface Si then 
u = iq and Eq. (1571) yields the required boundary conditions 

4>(Si) = 0(itj) = 0j with i = l,2. (40) 

The uniqueness theorems guarantees that such a solution is unique. Note that in this case the boundary conditions 0 1 and 
0 2 could have dimensions as well as the general solution 0(it). For instance, in electrostatics 0 (it) would have dimensions 
of potential. This strategy is quite similar to the use of Green’s functions to solve Poisson’s equation (or any other linear 
inhomogeneous equation), in the sense that Green’s functions only depends on the geometry as it is the case with the BHF’s. 

As a matter of example, we should keep in mind that the BHF’s given by Eqs. (15511551) for the concentric spheres are 
purely geometrical so far. If we want to solve a physical problem such as the electrostatic potential in the volume defined by 
b < r < a, when the two spherical surfaces are at potentials 0b and 0 a , such a potential is given by Eq. (1391) 

0 (r) = 0fa/i (r) + 0 a /2 (r) = 0 a + fi (r) [0 b — 0 a ] ; b<r <a (41) 


3 Formulas of capacitance in orthogonal curvilinear coordinates 

The BHF’s described in section [T] have an interesting application in determining coefficients of capacitance of certain geometrical 
configurations. Let us assume that the surfaces Si and Sj defined in section [j] are covered by electrostatic conductors and the 
volumen V is left empty. For example, let us assume for a moment that Sj contains to Si. We could have a conductor with a 
cavity such that the surface of the cavity is Sj and we could have another conductor inside the cavity whose external surface 
is Si as displayed in Fig. [T] Of course, it is possible to interchange the roles of Si and Sj and our developments are still the 
same. 
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Figure 1: Configuration of two conductors in which the inner conductor with external surface Si is enclosed within the cavity 
of an external conductor. The surface of the cavity is denoted by Sj, and the volume V is the region (in white) between both 
surfaces. We could interchange the role of the labels i j, and all our developments are still valid. 


We can determine the coefficients of capacitance for this geometry of capacitors by using the BHF’s in section [H It can be 
shown that the coefficients of capacitance for the previous configuration of conductors are given by jB] 


Gij — - 0 


V/j (iii) ■ n t dS 


where n; is a normal to the surface Si that points outwards with respect to the i — th conductor. We have 


(42) 


= ( 43 ) 

j 

with Qi, (f>i denoting the charge and potential on each conductor. Expression (l42jl is valid even if i = j. Nevertheless, we shall 
consider the case in which i ^ j henceforth. Since the functions fj take constant values on each surface, then V/j evaluated 
on each surface is perpendicular to such a surface. Moreover, it can be proved that [3] 


V/j {ui, v,w)-ni>0 for i ± j (44) 

such that V/j ( Ui,v,w ) is parallel (and not antiparallel) to in each point Ui,v,w on the surface Si. From these facts and 
using Eqs. (1721EB1) . as well as the fact that fj (u) only depends on u, we have 


V/j (m, v,w) ■ n 8 


V/j (' Ui,V,W ) • rq 


II V/j (ui,t>,w)|| 


1 dfj ( Ui ) 


1 dZ ( Uj ) 

h u {ui,v,w) du 


h u {ui,v,w) [Z (uj) - Z (m)\ du 


1 

h u (ui,v,w) [Z (uj) - Z (ui)\G(ui) 


(45) 


where we have taken into account that all functions involved in integral (1421) should be evaluated at a given point in Sj i.e. at 
u = Ui. Moreover, since dS is orthogonal to the unit vectors e„ and e w associated with the v,w coordinates, we have 


dS = \h v (ui, v, w) h w (ui,v,w)\ dvdw 

where we have taken into account that dS is positive definite. Utilizing Eqs. m the differential of surface becomes 

dS = \h u (ui, v, w) G (ui) H (v,w)\ dvdw (46) 

substituting (Tf5l) and (H51) in (H^l) the non-diagonal coefficients of capacitance yield finally 


r = _£o_ 

\Z ( Uj ) - Z ( Ui )\ 


\H [y,w)\dv dw ; with i,j = 1,2 and i^j 


Si 


(47) 


we emphasize that despite Eq. (H^|) is valid for any values of i and j, Eqs. (PT51) and (EPE1) are only valid for i ^ j. In order to 
obtain the diagonal coefficients we only have to use the properties [3] 

Gj -| T Gj 2 — 0 , Gj j — Cji Ca — Gij with i,j — 1,2 and i 4 J 
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or more explicitly 


Cn — C 22 — —C 12 — —C 21 


(48) 


It worths noting that expression (1471) is also invariant under the “gauge transformations” given by Eqs. (1251) . On the other 
hand, Eqs. (HTl) and (jlSl) leads to 0 0J 

Cn > 0 and C\j < 0 with i 7 ^ j (49) 

Finally, as a proof of consistency, we see from Eq. flgjl that the eigenvalues of the 2x2 matrix of capacitance are 0 and 2Cii. 
These results along with Eqs. (Hill) says that the matrix of capacitance is positive singular as it must be [5]. In terms of voltages 
Eqs. ([43l|48]) becomes 

Qi = C n (0i — 02) = C 22 (0i — 02) = ~Q 2 (50) 


Note that if S) refers to the internal conductor, Qi is the total charge on it. But if Si refers to the cavity of the external 
conductor that contains the other conductor, Qi is only the charge accumulated on the surface of the cavity. In addition, if 
the internal conductor has a cavity, the surface of such a cavity does not contribute in the surface integral to calculate the 
coefficients of capacitance [ 2 ]. 


3.1 Capacitance for the particular case of the scale factors 


As it was the case for the solutions of Laplace’s equation, the expression for the coefficients of capacitance is particularly simple 
when the system of coordinates obeys relation ECU) 


huh*. 


= 1 


(51) 


in that case, Eqs. (1271) provide simple solutions for Z (u) and H (v,w). Substituting (1271) in Eq. (j47[) . the latter becomes 


Cn = -1 


£ o 


dv dw 


7 L 3 


(52) 


or for the diagonal elements 


Cn = — Cij = - -- (f) dv dw ; with i,j = 1,2 and i ^ j (53) 

\Uj - Ui | J Si 

We saw that Eq. (ICTl) holds when we have a system of coordinates with cylindrical symmetry obtained with a conformal 
transformation from the cartesian coordinates. Nevertheless, we should recall that a cylindrical symmetry coming from a 
conformal transformation is only a sufficient condition. We only require the condition (1511) to guarantee the validity of Eq. 

(E3D- 

Observe that Eq. (l53l) has a structure similar to the capacitance of two parallel plates with “area” 

A e ff = (j) dv dw (54) 

located at a “distance” |Uj — Ui |. Note however that in generalized coordinates, the integral (1541) is not necessarily an area. For 
instance, if the two coordinates v and w are angular, this quantity is dimensionless. In the same way, the quantity \v.j — m\ is 
not always a distance. 


3.2 Capacitance of two concentric spheres 

Once again for the sake of illustration, we start calculating the well known capacitance of two concentric spheres with the 
values given in subsection 11.21 Thus, the radius of the smaller sphere is b and the radius of the cavity of the bigger sphere is 
a, so that a > b. Applying Eqs. (f3TT) and (l32l) in Eq. (l47l) and setting n = n = 6, rj = r 2 = a we get 


C 12 

C 12 


£0 

\ z O 2 ) - Z(n)\ 



| H (9 , ip) | dd dip = 


£0 ab 

C a~b) 



sin# dd 



4:TT£oab 

(a-b) 




| sin 0 1 dd dip 


it is straigthforward that C 21 = C 12 . Therefore the coefficients of capacitance yield the expected result 


Cii = — Cij = -; with i, j = 1,2 and * 0 j 

a — b 


( 55 ) 
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3.3 Capacitance of two concentric cylinders 

As for the case of BHF’s, it is more convenient to use the conformal polar cylindrical coordinates (p, ip, z) defined by Eqs. 
(1551) instead of the usual polar cylindrical system. Since for the (p, p, z ) system the condition (l5ll) is fullfilled the capacitance 
coefficients are given by Eq. (1551) 


Cu = 


£ o 


I Pj ~ Pi I Js, 

2-K£ 0 L 


dp dz = 


£ o 


I Pj Pi 
2n£ 0 L 


p2ir 

Jo ' 


ClJ \Pj-Pi\ | In (a/ 6 ) | 


dp 

hi = 


where we have assumed that the cylinders are located in the z interval [— L/2, 



1,2 and i^j (56) 

L/2] with L >> a. This is the usual result. 


4 Calculations of BHF’s and capacitances in more complex geometries 

As we said, the cases of two concentric spheres and two concentric cylinders were developed for the sake of illustration only. In 
this section, we shall calculate BHF’s and capacitance coefficients for more complex symmetries in which our method acquires 
all its power. 

4.1 BHF’s and capacitance on Prolate spheroidal surfaces 



Figure 2: Shell formed between two confocal prolate spheroids. 

Let us consider a shell formed between two confocal prolate spheroids as illustrated in Fig.[2j In order to find the BHF’s and 
capacitance coefficients, we use the Prolate spheroidal system of coordinates, which is obtained by revolutioning the elliptical 
bidimensional system of coordinates (£, rf) with respect to the major semiaxis of the ellipses. The most usual definition of the 
prolate spheroidal system of coordinates reads: 


x = a sinh £ sin 77 cos p ; y = a sinh £ sin r) sin p ; z = a cosh £ cos rj 

£ > 0, r?G[0 ,7r], a ~ I ^ G [°’ 27r ) 


where £, r) are the bidimensional elliptical coordinates, p is the azimuthal angle and d denotes the distance between the foci. 
The scale factors are given by 

he = h v = a\J sinh 2 £ + sin 2 rj ; h v = a sinh £ sin p (57) 

By means of the trigonometric and hyperbolic identities 


z 1 x 2 + y 2 

- 1 - - — 

a 2 cosh 2 £ a 2 sinh 2 £ 

z 2 x 2 + y 2 

a 2 cos 2 r] a 2 sin 2 77 


cos 2 77 + sin 2 77 = 1 
cosh 2 £ — sinh 2 f = 1 


(58) 

(59) 


it can be seen that the constant values of £ correspond to prolate spheroids, while surfaces of constant 77 are hyperboloids of 
revolution. Consequently, since we are interested in the configuration in which the surfaces are confocal prolate spheroids, the 
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coordinate that is kept constant is £. Then we shall assign (u,v,w) —> Combining Eqs. (1571) and (1771) we have 


hij h v 

h 


= hep = a sinh £ sin rj = G (£) H (r),(p) 


Such that a set of solutions for G (u) and H ( v , w) is given by 

G (£) = a sinh £ ; H (rj, ip) = sinrj 
and from Eqs. mm, the simplest solution for Z (u) yields 

ci£ 1 


(60) 


Z(0 = bn (tanh 


sinh£ 


di 


(61) 


substituting (RTT1) in (1771) we obtain 

/,(£) = 


Z(0-Z{& 


[zi^-ziti)] 


/,(£) = 


In (tanh — In (tanh l n ^ 


( tanh 4 


tanh - 


In 


(tanh y' j — In (tanh j n f 


tanh -j- 
tanh -it- 


i,j = 1,2 and i^j 


(62) 


where £, and £ 7 are the constant values of that coordinate on the surfaces S) and «Sj. Now for the capacitance we subtitute 


and (RTT1) in (l47l) to get 


C 12 — — 

Cl2 = — 


£o 


z^)-z(Zi)\r Sl 


-Y\^f s \ H (^^)l dl l d V 


/»27T pTl 

/ dip 

Jo Jo 


i In (tanh 4f ) - ± In (tanh 

the explicit calculation of C 21 gives C 21 = C 12 as it must be. So we finally obtain 

47ra£o 


sin 77 dy 


Ca — —Cij — 


In 


tanh -J- 


; i,j = 1,2 and i ± j 


(63) 


The reader can contrast the simplicity and brevity of the calculation in this subsection, with respect to the methods 
developed in the literature (see for instance Ref. [25] 1. Once the general formulas (1771 1771) and (HT1) are obtained, the 
application to this specific geometry is just a bit more than a direct replacement. This is also the case with other symmetries 
as we shall see. 

Finally, it is desirable to obtain the capacitance in terms of geometrical factors. Defining a,; as the major semiaxes of the 
prolate spheroid associated with £*, the capacitance (1671) becomes (see appendix [Cl) 


Cu — —Cij — 


47raeo 


In 


( ( a+ai)yjd 2 .-a 2 
\ (a+aj ) ■sja? —a 2 


; i,j = 1,2 and i ^ j 


(64) 


it is also shown in appendix [C] that the capacitance of two concentric spheres is obtained from the limit where a tends to zero. 

4.2 BHF’s and capacitance on oblate spheroidal surfaces 

Let us consider a shell formed between two confocal oblate spheroids as displayed in Fig. [3] Then, we should work with 
the oblate spheroidal coordinate system which is obtained by revolutioning the bidimensional elliptical coordinate system 
(£, rf around the minor semiaxis of the ellipses. The transformation of coordinates and scale factors yield 


x = a cosh £ cos rj cos ip ; y = a cosh £ cos y sin ip ; z = a sinh £ sin y 
he = h v = a\J sinh 2 £ + sin 2 rj ; = a cosh £ cos y 


£ > 0, |], a = 77 ; <P&{ 0,2tt) 


(65) 

( 66 ) 
(67) 
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Figure 3: A shell formed between two confocal oblate spheroidal surfaces. 


The trigonometric and hyperbolic identities 


a 2 skill £ a 2 cosiT £ 
x 2 + y 2 : 1 

a 2 cos 2 r] a 2 sin 2 77 


x 2 + y 2 2 , - 2 1 

- 2 — = cos 77 + sin 77 = 1 


= cosh 2 £ — sinh 2 £ = 1 


( 68 ) 

(69) 


show that the surfaces in which £ acquires constant values correspond to oblate spheroids, while the surfaces with 77 constant 
correspond to hyperboloids of revolution. Once again, £ is the variable that is kept constant in the surfaces we are working 
with. Equations (12^1) and (1M1) are combined to give 


hr,h 




= h v = a cosh £ cos 77 = G{£)H (r],<p) 
n t 

G(£) = a cosh £ ; H ( 77 , ip) = cos 77 

z (0 = ~ f \c d £=- 61111-1 (tanhg) 

a J cosh £ a 


Using Eqs. ([701 [71]) in (fTTl) we obtain 


fj (0 = 


Z (£) — Z (&) sin 1 (tanh£) — sin 1 (tanh£j) 


[Z (4j-) — Z (&)] sin 1 (tanh^j) — sin 1 (tanh&) 


I i,j = 1,2 and i 7 ^ j 


(70) 

(71) 


(72) 


and the capacitance coefficients are obtained from CDl EO) and dU) 

ry\^f s \ H dr i dL P = 


r . = ~ £ 0 

i^fe-)-^te)ir Sl 

c .. = _ c .. = _ 47r « g 0 _ 

I sin -1 (tanh £.,■) —sin -1 (tanh£j) 


£0 


4 sin 1 (tanh£j) — 4 sin 1 (tanh£j)| J 0 
■ ! i,j = 1,2 and i £ j 


p2iv p 7r/2 

dtp 1 

JO J-tt/2 


cos 77 dy 


(73) 


for instance, if we compare with the procedure in Ref. |26| , the advantages of our method are apparent. In terms of geometrical 
factors the coefficients Cj 3 read (see appendix [Cl) 


Ca — —Cij — — 


47raeo 


; i,j = 1,2 and i^j 


|cos £j — cos A ei| 

where ej, e* are the eccentricities of the ellipses associated with the surfaces S 3 and Si respectively. 


(74) 


4.3 BHF’s and capacitance of two confocal elliptical cylinders 



Figure 4: A shell formed between two confocal elliptical cylinders. The values 01 , 02 correspond to the major semiaxes, while 
61 , &2 denote the minor semiaxes. 
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This is an example in which the particular condition (EB1) is fulfilled. Let us consider a shell formed by two confocal elliptical 
cylinders as shown in Fig. [U The values b±, a±, correspond to the minor semiaxis and major semiaxis of the internal cylinder, 
while 62 , 02 are the semiaxes associated with the external cylinder. We assume that along the Z —axis the cylinders are located 
in the interval [-L/ 2, L/2] with L » ci 2 - We shall use the cylindrical elliptical system of coordinates which is obtained by 
projecting the elliptical bidimensional system on the Z— axis. The transformation of coordinates and scale factors yield 


x = a cosh £ cos 77 \ y = a sinh £ sin 77 


h% = h v = ay sinh £ + sin 2 77 ; h z = 1 
£ > 0 ; 0 < 77 < 27r ; —00 < z < 00 ; a = — 


The trigonometrical and hyperbolic identities 

x 2 


+ 


y 


a 2 cosh" £ a 2 sinh 2 £ 

2 2 

x y z 


a 2 cos 2 77 


9 • 2 
a* sm 77 


cos 2 77 + sin 2 77 = 1 
cosh 2 £ — sinh 2 £ = 1 


(75) 

(76) 

(77) 

(78) 

(79) 


say that surfaces with constant £ are confocal elliptical cylinders and surfaces with 77 constant are confocal hyperbolic cylinders. 
Thus £ is the coordinate that takes constant values in the surfaces we are considering. Equation (l76l) shows that in this case, 
the condition (l26l) is satisfied. Consequently, the BHF’s are given by Eq. (l28l) that in this case becomes 

fj (£) = J—% ; bi = !, 2 and i ^ 3 (80) 

si S i 

and the coefficients of capacitance are given by Eq. (15^1) that in our case gives 


II 

<5 

—C 12 — 

Ca = 

P 

11 


£0 


I £2 - £i| Ts 1 
2 tt£ 0 L 


di 7 dz = 


£0 


16 - 6 


16 - £i| Jo 

i,j = 1,2 and i^j 


/»27t /*L/2 

/ dp dz 
Jo J —L/2 


(81) 


Note that the integral has the same value if we integrate over S 2 , it guarantees that C 12 = C 21 . This example shows that 
under the condition (1^1) the results are much simpler. In terms of geometrical factors Eq. (l 8 lT) becomes (see appendix [Cll 


Cu — —Cij — 


2tt£ 0 L 


In 


i + -\J af—a 


aj+y/af-a 2 

the limit a —> 0 , reduces to the case of two concentric cylinders. 


hj = 1,2 and i ^ j 


(82) 


4.4 BHF’s and capacitance of two eccentric cylinders 



Figure 5: Two eccentric cylinders C\ and C 2 with radii rq and r 2 . The cylinder C\ is contained in the cylinder C 2 , and their 
centers are separated by a distance D. 

Let us consider a shell like the one in Fig. [5] The external cylinder has a radius rq, and the internal one a radius r\. The 
centers are located at the points (x ci ,0) and (a; C 2 ,0), and are separated by a distance D. The cylindrical bipolar coordinate 
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Figure 6: Bipolar coordinate system (in two dimensions) with foci F\ = (0, — a) and F 2 = (0, a). It is shown that curves 
associated with positive (negative) constant values of r are circles that contain the focus F 2 (Fi). When |ti| > |t 2 | the circle 
G 2 contains the circle C\. 


system is the appropriate for this geometry. The (bidimensional) bipolar coordinate system posseses two foci F\ and F 2 , 
usually located at (—a, 0) and (a, 0) respectively. The transformations from cartesian coordinates and scale factors are given 
by 


sinh r 


x = a 


h r , = h T = 


cosh r — cos a 
a 


; V = a 


sin a 


cosh t — cos a 
h, = 1 


z = 21 


cosh r — cos cr 
0 < a <2i r ; —00 < r < +00 ; —00 < z < +00 

The distance from each focus to an arbitrary point (x, y ) in the plane yields 

d\ = (x + a) 2 +y 2 ; d 2 = {x — a) 2 + y 2 

Thus, a represents the angle between the lines that join each focus with a given point (x, y) in the plane, while r is given by 


r = In 


the curves (in the plane) with a constant correspond to non-concentric circles that intersect in the two foci (see Fig. 0. The 
curves with r constant correspond to non-intersecting circles and with different radii (see Fig. [6]). Two curves C\ and C 2 with 
constant positive values n > T 2 correspond to two non-concentric circles such that both circles contain the focus F% and the 
circle C 2 contains the circle C\. Adding the z coordinate, we obtain the desire configuration of two non-concentric cylinders!!. 
Therefore, we assign ( u,v,w ) —> ( t,<j,z ). Note that this system of coordinates satisfy condition C()l) so that the solutions for 
the coefficients of capacitance and BHF’s are given by Eqs. (E51) and (15^1) 


fj ( T ) = 
C u = 


T ~ Ti 


£ 0 


; i,j = 1,2 and i ± j 


dcr dz = 


£ o 


Ca — ~Cij — 


t 3 ~ T i\ JSi 

2tt£qL 


Ti - Ti 


/*27T 

■ / da 

Jo 


f-L/2 


dz 


-L/2 


; hj = 1)2 and i ^ j 


(83) 


(84) 


\Tj ~Ti | 

Substituting Eq. (11221) of appendix [C] in Eq. (1M|) , we obtain the coefficients of capacitance in terms of geometrical factors 

2tt£qL 


C 1 — —C 1 — 


when D — > 0 we recover the expression for two concentric cylinders. 


(85) 


$We can choose negative values of r such that the circles contain the focus F \. 
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5 Conclusions 

We solve Laplace’s equation for some geometrical configurations in which there exist a system of coordinates (u, v, w) that can 
be adjusted to the geometry. Specifically, we assume that the boundary consists of two closed surfaces Si and S 2 so that one 
of the surfaces is included into the other, and in which the coordinate u takes constant values u\ and U 2 on each surface S± and 
S 2 . We solve Laplace’s equation within the volume V between both surfaces. Besides, we take specific boundary conditions 
in which the solutions take values that are either zero or unity on the surfaces. The solutions of Laplace’s equation with such 
specific boundary conditions are called Basic Harmonic Functions (BHF’s). It worths remarking that the BHF’s only depend 
on the geometry, and their solutions via separation of variables are particularly simple. The expressions for the BHF’s are 
given in terms of generalized curvilinear coordinates in such a way that we can sweep a variety of systems of coordinates to 
obtain the BHF’s associated with the given symmetry. 

Further, we exploit the purely geometrical BHF’s in two ways. On one hand, we can determine a more general class of 
Harmonic Functions (HF’s i.e. solutions of Laplace’s equation) for the same geometry of the BHF’s, but in which the functions 
take arbitrary constant values on the boundaries. On the other hand, we calculate the coefficients of capacitance when the 
surfaces that form the volume are covered with electrostatic conductors. Expressions for either HF’s and for coefficients of 
capacitance are given again in terms of generalized curvilinear coordinates, such that we can obtain these quantities for multiple 
symmetries by simple replacements. Several examples with different symmetries are provided. 

In addition, appendices lAl and iBl show natural extensions coming from our present formulation. In appendix 0 expressions 
for the electric field between two conductors and the charge densities on each of their surfaces are written in generalized 
curvilinear coordinates, while in appendix [Bj formulas for the capacitance coefficients are extended to the case of conductors 
in sucessive embedding. 

It worths emphasizing that after obtaining the general formulas for BHF’s, capacitance coefficients, electric fields and 
surface charge densities in terms of generalized orthogonal curvilinear coordinates, calculations of these observables for each 
symmetry are straightforward. After checking Eqs. (1551 15511771) and Eqs. (155115U1) we can see that all those observables are 
obtained by simple combinations of the functions h Ul G ( u ), H (v. w) and Z ( u ). The first three are basically input parameters, 
and the latter is obtained from 1 /G (it) by a direct integration. Comparison with the calculations of the same quantities with 
other methods given in the literature, shows the simplicity and brevity of our method as it was emphasized in several of our 
examples. 

Note that despite with the traditional strategies Laplace’s Equation can be solved in many cases via separation of variables, 
in most of such cases a series expansion in certain special functions is usually required in order to adjust the boundary 
conditions. By contrast, the BHF’s obtained via the formulas (15511551) lead in many problems to closed solutions. Therefore, 
solutions of Laplace’s equation with equipotential surfaces can often be written (with our method) in terms of closed well- 
defined functions. Though solutions in series are exact in principle, when we handle them in practice we should truncate the 
series somewhere, leading to a considerable expansion of numerical errors, especially if the series converges slowly (sometimes 
problems of convergence also arise). 

On the other hand, it is well known that solutions of Laplace’s equation can be treated with the geometrical formalism of 
Green’s functions. Nevertheless, it could be noticed that finding BHF’s (which are also purely geometrical) is in general much 
easier than the calculation of the Green’s function on the same geometry. Finally, it is possible to extrapolate the method 
developed here to other linear homogeneous differential equations. In the same way that we can associate a Green’s function 
with each linear differential operator, we can define BHF’s for each linear differential operator. 
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A Electric field and surface density in a configuration of two conductors 


Once again we use i,j = 1,2 and i ^ j throughout this section. If <j>i and <pj are the (constant) potentials on conductors i and 
j respectively, the potential in the volume V between the surfaces S, and Sj is given by Eq. (1551) 


f — ^i/it 0j fj 

on the other hand, from the property _/) + fj = 1 we obtain 

V/i = -V/i 


( 86 ) 

(87) 


combining Eqs. (15511571) and the fact that the BHF’s only depend on u, we can find the electric field within the volume V in 
terms of curvilinear coordinates 


E = -V<t> 


-oX.f, - faVfa = {fa - fa) S7fj 


fa ~ fa dfj (it) 
h u {u,v,w) du 


(u,v,w) 
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where e u (u, v, w) is the local unitary vector associated with the u coordinate. Using Eqs. 

fa — 4>j dZ (u ) 


we have 


E (u,v,w) = 


-e u (u,v,w ) 


h u (u , v, w) [Z (Uj) — Z (itj)] du 
(f>i — cj)j 

E (u, v, w) = — - - —-- e u (u,v,w) ; i,j = 1,2 and * ^ j 

h u (u,v,w) [Z (uj) - Z (m)\ G(u) 

On the other hand, the surface charge density induced on the surface Si of a perfect conductor is given by [2] 
c7 (Si) = a (m, v, w) = £oV <j> ( Ui , v, w) • n* = e 0 (fa - fa) V/j (it,, u, w) ■ n, 
where we have used Eqs. (1551) and (1571) . Finally, by using Eq. (H5l) we have 

£o 


( 88 ) 


a(ui,v,w) = 


(fa - fa) ; i,j = 1,2 and i ^ j 


(89) 


h u (ui,v,w) [Z (ufa - Z (ufaG (u^ J 

As a proof of consistency, let us obtain the total charge on the surface Si. By using Eqs. (15511751) and (14711751) we obtain 

dS 


Qi = (b <j(ui,v,w) dS = e 0 (fa - fa) <h 

JSi JSi 

= E 0 (fa-fa)(j) 


h u (ui, v, w) [Z (uj) - Z (ui)]G(ui) 


h u (ui, v, w) G(m) H(v,w) 


dv dw 


— (fa fa)' 


I Z(uj) - Z (m )| fs. 
Qi — Ga ((/>i 4*j) 


h u (ui, v, w) [Z (Uj) - Z (ui)]G(ui) 

(j) \H (v,w)\ dv dw = —Cij (fa — fa) 


(90) 


which reproduces the correct expression (l50l) for the total charge on Si. If we compare expressions (1551) and (1901) . the term 
in the absolute value of Eq. (1551) behaves like an “effective capacitance per unit area”. Effectively, such a term is purely 
geometrical, as is the term Cu. 

A.l Electric fields and charge densities for two concentric spheres 

For the two concentric spheres of Sec. 11.21 we substitute Eqs. (125115111521) in Eq. (1551) to obtain 

(4*b 4*a) 


E (r) = 


E = ab(fa fa) ^ 
(a — b) r z 

and substituting (155115T11551) in Eq. (1551) we find 


n) = 


4*j ) 

(fa - 

fa) e r 

hr [Z fa) - Z (n)} G (r) 

1 _1_ 

r 2 ri 

far 2 ) 


far 2 ) [i - i] (-r 2 ) 


(91) 


£o 


hr (r%) 


fafa) 


(fa - fa) = 


Tj Eq ((f)i (fa ) 


Ti r f - r,- 


/ \ h Eq faa fa) /., El Eq (fa (j) a ) 

G W = 7—77-U— ; cr (o) = T- 


(92) 


a (a — b) ’ b (a — b) 

integrating these densities on S a and Sb we obtain the total charge on such surfaces. Further, taking into account Eqs. m 
EU) we obtain the consistency of all these results 

47re 0 a& (fa - fa) 


Qa — Qb — 


= On (fa -fa) ] E (r) = 


Qb 


(a — b) ’ 47T£ 0 r 2 

A.2 Electric fields and charge densities for a prolate spheroidal shell 

For a prolate spheroidal shell, we obtain the electric field within the shell and the charge density on the ellipsoidal surfaces by 
substituting Eqs. (l57l l60l) and (15T1) in Eqs. 


E(£,7?,¥>) = 


G (Ui, Tj, fa) = 


Hi faj 


a fa sinh 2 £ + sin 2 ry In ^tanh faj — In ^tanh 


sinh£ 


- e « (€,V,<P) 


£ o 


a fa sinh 2 fa + sin 2 In ^tanh faj — In ^tanh faj 


sinh fa 


(fa -fa) ; i,j = 1,2 and i fa j 
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the magnitude of the electric held and the surface charge densities are independent of the azimutahl angle <p as the symmetry 
suggests. Calculation of electric fields and surface charge densities in other complex geometries is also straightforward. 


B Chains of embedded conductors 



has only an inner part (the surface of the cavity). Conductors 2 and 3 have an inner part S\ b ' 1 and an outer part S\ 


:(“) 


Let us assume that we have a set of N + 1 conductors succesively embedded as indicated in Fig. [7] We label them 
as 1, 2,..., N + 1 from the inner to the outer. We can see that for k = 2,..., N the surface Sk of the k — th conductor has an 
inner part that we denote as s [ b ^, and an outer part that we denote as s!, a ' > . However, for Si we only have outer parl{§, while 
for Sn+i we only have an inner part. Nevertheless, we still preserve the inner and outer notation for Si and Sn +i by writing 
Si = S[ a ' > and Sn+i = S$ +1 . The volume 14 with k = 1,..., N is defined between S ^ and S^_ v 

A set of conductors in succesive embedding has some particular properties (see appendix C in Ref 0). First, the only 
non-zero coefficients of capacitance are the diagonal ones Ca and the terms of the form Ci t i± i. Of course, for i = 1 there is 
no term of the form Cj^—i, and for i = N + 1 there is no term of the form Ci ^+The diagonal terms are related with the 
non-diagonal ones as follows 


Cu — — (C^j-i + Cj^+i) ; i — 2,...,N 

On = —O12 ; Ojv+i.jv+i = —C^jv+i (93) 

now, owing to the symmetry of the matrix and the properties ( 1931 ) . it is enough to calculate the non-diagonal terms of the 
form Cj.j+i for i = 1,..., N. It can also be proved that jj] 


Ci li+ i = -e 0 (f V/ i+1 (u| a) ) • n ( “ ] dS ; 

J S) 




(94) 


that is, only the outer part of the surface Si contributes to the surface integral in the calculation of CVj+i . It worths pointing 
out that the properties described above are valid for any geometry of the conductors, the only condition is that they should 
be in succesive embedding. 

Now let us add the hypotheses we have been working with. That is, that we have a curvilinear coordinate system ( u , v, w) 
such that the coordinate u acquires a constant value Up X ^ on each surface S^ with p = 1,2, ..., N + 1 and A = a, b. 

First of all, we want to calculate the BHF’s within each volume V). By definition of the BHF’s, they satisfy the boundary 
conditions 

fi (X (A) ) = 1, fi (sj A) ) =0 ; ijtj and A = a,b (95) 

It could be shown that within the volumen V) the BHF’s have the following properties 

fi (Vi) + f i+ i (Vi) = 1 ; fi+k+i (Vi) = fi-k {Vi) = 0 ; i = l,2,...,N , k > 1 


§Let us recall that even if the most interior conductor has a cavity, we can exclude its surface for calculations. 
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therefore, it is sufficient to calculate /*+1 within the volume V). Our boundary conditions for the BHF f t +\ within Vi are given 
by 


f i+1 (s[ a) ) = 0 ; / i+ i(sj5i) 


= 1 


with the same analysis done in section|l] we see that the functions fi+\ only depends on u, such that these boundary conditions 
can be written aqj 

' = 1 (96) 


fi +1 (*4 o) ) = ° ; fi+i («i+i) 


now, taking into account that in this context the volume Vi in which we are evaluating the BHF’s is delimited by the two 
surfaces Sj a ^ and we can use the same arguments that led to Eqs. (l22l) with the assignments 


.(«) 

x i > 


r (b) 

l i+1 


i —t i , j —> i + 1, Ui 

From these facts the function fi + \ evaluated at any point within the volume V); yields 

Z (u) - Z (u^) 


fi +1 (V) = 

With a similar procedure we can show that 

/*+1 ( Vi + 1 ) = 


(«&) - * (tt, (B) ) ' 

Z(u)-Z(u% 2 ) 
z (uj+i) - Z (u[%) 




(97) 


(98) 


; i = 1,..., N — 1 


f l+1 (V l+1+k ) = f i+1 (Vi_ k ) = 0 ; k > 1 


(99) 


( 100 ) 


On the other hand, the coefficients of capacitance can be obtained from (1M1) and [Ml with the same arguments that led from 
QID to (03). Thus taking into account the assignments (1971) and the fact that only the outer part of the surface Si contributes 
to the surface integral in the calculation of we obtain 


Ci,i+i — — 


/ 777 N°° - / 4, J H (v,w)\dvdw ; i = 1,2, ...,7V 


( 101 ) 


B.l Succesive embedding of concentric spheres 

A simple example consists of examining the succesive embedding of N + 1 concentric spherical conductors. In that case u = r 
and we define a set of radii in the form 


a 3 <b 2 < a 2 <b 3 < a 3 < ... < < &jv < 


( 102 ) 


note that bi refers to the radius of the cavity of the i — th conductor while ai refers to the radius of the external surface of the 
i — th conductor. Hence, ai — bi is the thick of the i — th conductor. 

Using (1M1) and (1571) . we obtain the solution for /*+1 in a volume V) (i.e. between the radii a,; and bi + 1 ) 


f (] = Z ^ Z ( r * 0 _ Z(r)-Z (a^ _ y - jr 

Z (r&) - Z (r<“>) z ( & Wi) - Z (a,) ^ ± 

fi+ 1 (r) = for ai<r < b i+1 

ai - b l+ i \r ai) 

and substituting Eqs. (1321 l3ll) in Eq. (I101D the coefficient of capacitance C)y +\ yields 


(103) 


Ci.i+i 


Ci,i+1 


£o$g(») |77 (#,</?) |dd dtp 

\ Z ( r i+i) — Z ( r | 0) ) 

47T£ 0 aA+i 

^j+i ai 


£o sin $ d9 dp 
\Z (b t+1 ) - Z (a x )\ 


47T£ 0 



(104) 


'Note in passing that the fact that only fi and fi+i are non-zero in Vi could be understood with the boundary conditions on Vi. For instance, 
the boundary conditions on Vi for fi—i would be fi—i = 0 and fi—i = 0, by uniqueness the only solution within Vi is fi—i (Vi) = 0. 
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C Curvilinear coordinates and geometrical factors 

We have obtained coefficients of capacitance in terms of the generalized curvilinear coordinates ( u,v,w ) that are adapted to 
their symmetries. Notwithstanding, since such coefficients are purely geometrical in nature, it is more useful to obtain formulas 
for them in terms of geometrical factors (radii, major semiaxes, focal distances, etc.) instead of generalized coordinates (it, v, w). 
Hence, in this appendix we rewrite the curvilinear functions obtained for the capacitances in terms of geometrical factors. 


C.l Prolate confocal spheroidal shell 

Recalling the equation of an ellipsoid of revolution (revolved around the Z— axis) in cartesian coordinates, we have 

z 2 x 2 + y 2 

9 -T 79 


(105) 


for a prolate spheroid, a* refers to the major semiaxis and bi to the minor semiaxis of the ellipsoid of revolution generated by 
Comparing equations (1551 11051) and taking into account that cosh£ and cosh£ are positive for £ > 0 we have 


Gtj Oj 

cosh £ = — and sinh£ = — 
a a 


we shall also take into account that 



by using a hyperbolic identity, and applying Eqs. (110611 we have 


tanh — 
2 

tanh — 
2 


sinh & _ bj/a _ b t 

1 + cosh & 1 + ^ a+ a,i 

\/ a i - « 2 

a + ai 


(106) 

(107) 


(108) 


note that the focal distance d = 2a is the same for all ellipses by definition. Therefore, it is convenient to put the expressions 
in terms of the major semiaxes a,;, aj of each ellipsoid and the (common) focal semi-distance a. Substituting (110811 in Eq. (16811 
we obtain the coefficients of capacitance for a prolate confocal spheroidal shell 


Cu — -Cij — 


47raeo 


In 


(a+ai)y/a 2 -a 2 
(a+aj)^ a 2 —a 2 


(109) 


it is interesting to check that despite equation (11091) is not defined at a = 0, the limit a —> 0 exists and reduces to the case of 
concentric spheres as it must be. Such a limit can be found by using L’Hopital’s rule. Assuming a 3 > at we have 


lim 

CL —^0 


47T£o 

Cu 


lim 

CL — >-0 


47T£ 0 



In 


lim ■ 

a —>0 


lim 

CL — >-0 


(a + ai) Ji 


a) - a -> 


— In 


(a + a,j ) yja 2 — a 2 


n 2 n 2 (a+a»)a 

a i sM^ 


J a 2 _ a 2 __ (a+ajW 


lim 

a—>0 

{a 3 - ai) 


(a + tn) yja 2 - a 2 {a + a 3 ) yj a 2 - a 2 

(ai — a) (a 2 — 2a 2 — a^a) — (a 3 — a) (a 2 — 2d 2 — a 3 a ) 


(a 2 - a 2 ) (a 2 - a 2 ) 




by comparing with (15511 we see that the limit of concentric spheres is obtained appropriately. 


C.2 Oblate confocal spheroidal shell 

In the case of an oblate ellipsoid of revolution (revolved around the Z— axis), its equation in cartesian coordinates, yield 

z 2 x 2 + y 2 

6? «? 


( 110 ) 
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comparing Eqs. (11101 1551) and using Eq. (11071) we have 


a 2 sinh 2 & 
tanh£i 


b 2 ; a 2 cosh 2 & = a 2 => tanh £* = — 

di 



V a i - « 2 


di 


( 111 ) 


where we have taken into account that a/cq = e*, where denotes the eccentricity of the ellipse associated with the coordinate 
Substituting (11111) in (1751) we have 


C u = 


Anaeo 


sin 1 (\/ 1_e ?) “ sin 1 (V 1 - e *) 

using the identity cos -1 2 = sin -1 \/l — x 2 for 0 < x < 1 we get 

47raeo 


C u = 


cos -1 e, — cos _i e 


-i, 


which coincides with (1751) . Once again, the limit with a —> 0 reduces correctly to the case of concentric spheres. 

C.3 Two confocal elliptical cylinders 

By assuming that the major semiaxis at goes along X and the minor semiaxes goes along Y we have 


2 2 
x y z 


b 2 


= 1 


( 112 ) 


and comparing with (1781) , we obtain 

we shall use the identity 


tanh£i = — 

di 


tanh 1 x = — In ^ + 


\ — x 


by applying (11141) and (1 10 71) in Eq. (11131) yields 


; |*| < i 


(113) 

(114) 




so that 




(a,i + bi ) 2 


5- 


{ai + bi) 


= In 


di + bi 


&_fj=ln| ^Ty l = 1 ” 


k + \J d 2 


di + a/ a 2 — a 2 


(115) 


substituting (11151) in (15TT) we obtain Eq. (1571) . This result is clearly independent on the assumption that the major semiaxes 
goes along X. 


C.4 Two eccentric cylinders 

We shall find the coordinate quantity T 2 — t\ in terms of geometrical factors: that is in terms of the radii rq and ?~2 of the 
cylinders and the separation D between their centers. We shall work in the two dimensional space with coordinates (r, a). We 
shall start from the identities 


x 2 + (y — a cot a) 2 
y 2 + (x — a coth t) 2 


. 2 

sin a 


sinh 2 r 


( 116 ) 
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and the equation in cartesian coordinates for a circle of radius r centered at ( x c , 0) 



(x - x c ) 2 +y 2 =r 2 

(117) 

coniDarine (|116p with (I117|) it 

can be shown that 



a = r sinh r ; x c = r cosh r 

(118) 

from Eqs. (Ill8l) we have 


a = r\ sinh t\ = r 2 sinh T 2 

(119) 


D = x C2 — x Cl = r 2 cosh T 2 — r\ cosh t± 

(120) 

squaring Eq. (I120D we have 


D 2 = r\ cosn 72 + r\ cosli T\ — 2r\r<i coshri coshr 2 

(121) 


and using the hyperbolic identity cosh 2 r — sinh 2 r = 1, as well as Eq. (11191) . then Eq. (I121D becomes 

D 2 = r\ + r\ + r 2 sinh 2 T 2 + r\ sinh 2 n — 2 r\T 2 cosh t\ cosh T 2 = r\ + r 2 + 2 (ri sinh t \) (j "2 sinh 72 ) — 2 rir 2 coshn cosh T 2 

= 7'| + r 2 + 2 rir 2 (sinh T 2 sinh n — cosh n cosh T 2 ) 

D 2 = ^ + r 2 — 2nr 2 cosh (r 2 — n) 

from which we obtain the desired relation 

T 2 - n = cosh" 1 ^ + 2 ^ r ~ ^ ( 122 ) 
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